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Abstract 

In a semiconductor-based double quantum well (QW) coupled to a degree of freedom with an internal 
dynamics, we demonstrate that the electronic motion is controllable within femtoseconds by applying ap- 
propriately shaped electromagnetic pulses. In particular, we consider a pulse-driven A^Ga^As based 
symmetric double QW coupled to uniformly distributed or localized vibrational modes and present ana- 
lytical results for the lowest two levels. These predictions are assessed and generalized by full-fledged 
numerical simulations showing that localization and time- stabilization of the driven electron dynamics is 
indeed possible under the conditions identified here, even with a simultaneous excitations of vibrational 
modes. 
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The impressive advance in fabrication, characterization and control of nano-scale systems [QQ| 
has fueled wide-ranged studies on the fundamental behavior of quantum systems and their utiliza- 
tion for practical applications, e.g. as elements in nano-circuits [2], nano-oscillators (3]|, nano- 
magnets [4] and opto-electronic devices 01. A similar rapid progress was also achieved in optics. 
Light sources are currently available in a wide range of intensity, duration and pulse shapes flU, 
offering new possibilities for the temporal and local manipulations of electronic properties of 
nanostructures. A paradigm example has been the ultrafast control of quantum states in a double- 
well (DW) structure which is of relevance for quantum computing application @. In this respect, 
three aspects are important: i) the swift electron localization in one of the wells starting from an 
arbitrary initial state, ii) maintaining this localization for a desired time, and iii) switching of the 
localization between the two wells. From a practical point of view, it is also essential to address the 
role of the coupling to the environment. In this work we study the electron dynamics in a typical 
Al^-Ga^As-based QW coupled to vibrational modes with the aim of controlling and sustaining 
the electron motion on a femtosecond time scale via electromagnetic pulses. The DW potential 
barrier V con f(x) with a well separation x confines the conduction electrons in the x-direction 
and is well modelled by V con f(x) = 16 (Vb/xq) (x 2 — Xq/4) 2 . Typical values for this system are 
Vb ~ 240 meV and x « 10 nm. The effective mass is assumed constant across the structure 
m m 0.066m e (as for GaAs, and m e is the free electron mass). The tunnel splitting hVt = E 2 — E\ 
where Ei and E 2 are the two lowest eigenenergies, is smaller than the separation from the third 
excited state, hence a consideration within a two level system (TLS) approximation is useful for 
low excitations and will be employed below for the analysis of qualitative trends. In addition we 
couple the electron dynamics to a degree of freedom (y) with an internal vibrational dynamics 
with the Hamiltonian Hr = p 2 /2m + ky 2 /2 and a coupling potential He- The full Hamiltonian 
reads then H = H + He + Hr. Here H is the bare electron term; the coupling has the form 
He — Vo{y/yo) v d x ) with the function v c (x) describing the interaction strength along the struc- 
ture, y stands for the typical length scale associated with the oscillator potential. This problem 
resembles the rotational dynamics coupled to vibrations [7] and also inelastic scanning tunneling 
experiments [8] . We will discuss two cases: (i) a uniform harmonic force acting on the electron, i.e. 
the coupling potential is v c (x) = x/xq, and (ii) a localized coupling v c (x) = exp[— (x — x c ) 2 /2\ 2 ]. 
Physically such vibrational (phononic) coupling may arise due to inelastic scattering of the tun- 
neling electron in the DW, e.g. from pinned impurities that then vibrate around the minimum of 
the pinning potential. 
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A particularly efficient and swift mean for steering the electron dynamics in DW is the use of 
specially shaped femtosecond light pulses such as the so-called half cycle pulses (HCPs) ll9l[T0l. 
These are strongly asymmetric pulses with one short and strong half-cycle followed by a very 
weak and long half-cycle of opposite polarity. A key point in this control scheme is that the pulse 
duration has to be much shorter than the characteristic time of DW system (T = 2n/VL) [TTOl . In 
this case, within the dipole approximation, the action of the pulses can be viewed effectively as 
impulsive kicks [11] delivered at times t k with a strength Ap fc = J E k dt where E fc is the electric 
field amplitude of the pulse and t k is the time when E fc reaches its peak, meaning that such a pulse 
sequence acts in effect as E fc = Ap k 5(t — t k ) and couples to the electron via H L (t) = xE(t) lfl~2ll . 
This type of pulses can localize the electron within femtoseconds whereas the localization time 
for harmonic pulses is about few picoseconds lfl~3l . How effective this scheme would remain if the 
tunneling electrons are coupled to an external degree of freedom is to be clarified here. Physically 
one would expect a similarly effective control scheme when the electron dynamics is coupled to 
an external bath as long as the time scales involved are longer than the pulse duration. However, 
the whole process is a coupled dynamical one with a non-trivial dependence on the properties 
of the external degree of freedom. We start thus from the two-level approximation. The case 
of zero vibrational coupling has been analyzed analytically lfT0l[l4l . The electron in the ground 
state is completely delocalized across the heterostructure. In what follows we assume that in 
the ground state the electrons are decoupled from the vibrations. The vibrational coupling acts 
on the excited states, i.e. as soon as the external pulses are applied. Starting from this state, 
in order to localize the electron within the left well, it is sufficient to apply a single 7r/2-pulse 
in x-direction, i.e. E{t) = Ap 5(t) with a = 2fiAp/h = ir/2, where /i being the transition 
dipole between the first two states and a can be viewed as the "pulse angle". After acquiring an 
additional momentum shift Ap the electron tunnels to the left well within the time tl = To/4. As 
the system undergoes free Rabi-oscillations, the probability of finding the electron in the left well 
Pl is given by Tz(i) = cos 2 (fit/2) . Within the half-period of the Rabi oscillation the electron 
is then mainly localized in the right quantum well. To maintain the localization, one may kick 
periodically the electron back by an appropriate laser pulse. Applying the 7r-pulse (a = n) to the 
free-evolving system at t — to brings the system back to the localized state at t = 2t so that 
-Pl(O) = P L (2t ) = 1. Therefore, the system has the periodicity T of the array of laser pulses 
with T = 2t . It is clear that the averaged localization will be higher if the interval T between two 
consecutive pulses is shorter. 
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To proceed further, we map the electron Hamiltonian H onto the spin-like TLS Hamiltonian 
Htls — —hVLa z /2 (with o~j being the j Pauli matrix). The mapping can be extended to incorpo- 
rate the vibrational coupling by substituting the position operator with —x a x /2 (with eigenstates 
localized in the left, right potential minima). Applying this rule to the coupling term v c (x) leads 
to the so called spin-boson-model (see [15J and ref. therein). According to the theory of small 
polarons lfT6l . the remaining vibrational degree of freedom y is incorporated as a renormaliza- 
tion of the tunneling amplitude Q — >• f2 = M n ^ n Vt. Thus, the polaron transformation decouples 
the electron from the vibrational degrees of freedom, so that the remaining renormalized elec- 
tron Hamiltonian can be handled analytically. We will use this model to make predictions for the 
modified control mechanism and compare with the full numerical exact result. 

(i) Uniform coupling v c (x) = x/x . The effective electron Hamiltonian reduces to the simple 
form Htls — —hVta z /2 with the renormalized tunneling amplitude Vt. For the ground state vi- 
brations, the renormalization is given by the Franck-Condon factor M = exp[— V 2 / (2h 2 co 2 h )}. 
Thus, the effect of the vibrations is to slow down the electron dynamics. All control strate- 
gies remain unchanged, but the time for achieving an initial localization increases to T£ = 
(7r/2f2) exp[VQ / (2h 2 uip h )]. For an application of this model to the DW coupled to localized 
phonons, we may choose the Debye temperature as the maximum phonon energy, i.e. huj vh = k B Q 
IfTTl [T8ll , and measure the coupling constant V in units of hui p h- In the phonon-free system, 
the characteristic time is T ~ 80 fs. Due to electron-phonon (e-ph) coupling, it increases to 
T ps 110 fs. 

(ii) Localized coupling v c (x) = exp[— (x — x c ) 2 /2\ 2 }. Substituting again x — > —x a x /2, and 
evaluating the matrix exponent yields 

x 2 /4 + x 2 c ~ 



y 

H c -> V — exp 

yo 



(cosh£ — a x sinh£), (1) 



2A 2 

with £ = xox c /2\ 2 . If the coupling to the phonon degree of freedom is centered around the 
potential barrier, i.e. x c = and therefore sinh £ = 0, He will no longer depend on the electron 
term. The two degrees of freedom are then decoupled. For the full spectrum system, this means 
that the phonon effect plays a minor role because the wave function nearly vanishes at x = 0. 

The worst case scenario in terms of invalidating previous localization schemes is the e-ph in- 
teraction acting in one of the quantum wells, since the resulting problem has now the maximal 
asymmetry. For the TLS model (after the polaron transformation), this asymmetry leads to an 
additional term H TLS = —Mla z /2 =f Kk<j x /2 for x c = ^fx /2. The parameter k is quadratic in 



V : k = q(X)VQ /h 2 oj p h, where q(X) = 1 — e~ x o/ 2A2 is a prefactor describing the spatial extension 
of the e-ph interaction along the structure. For this situation we inferred the following results: 

(a) Localizing the electron in the left well: Being no longer an eigenstate of the uncoupled Hamil- 
tonian, the symmetric ground state of the DW potential oscillates and therefore gains the momen- 
tum T$Pi f° r %c = T%o/2 an d t < T /4. The corresponding pulse angle a must be modified 
according to a = n/2 =)= Afix with Afli = 2fx5pi/h. 

(b) Maintaining the localization: Due to the asymmetry, the electron tunnels back from the left 
well with an increased (a decreased) momentum p ± 5p 2 for x c = xq/2 (x c = —Xq/2). Note that 
this additional momentum 5p 2 depends on the time. In order to incorporate the dynamics arising 
from 5p2, the strategy is now to use the same pulse sequence, but to adjust the pulse angle to 
a = 7T =f A/3 2 with Af3 2 (t) = 2/i5p 2 (t)/h. 

The angles A/3i and A/3 2 , or, equivalently the transferred momenta, can be calculated using the 
TLS model. 

Alternatively, one can also compensate the asymmetric vibrational coupling by a constant electric 
field E oriented in the x-direction. In the presence of the electric field = ^fhn/2fi, the TLS 
term proportional to a x is exactly canceled. Thus, if the value Eq^ is tuned appropriately, the 
control strategies remain the same as in the case of zero coupling. 



In order to compare the optimal values obtained from the TLS model with the exact re- 
sults, we performed a full numerical propagation based on the two-dimensional split-operator 
approach [fT9l . For definiteness we consider x c = —x /2, the case x c = x /2 can be treated 
analogously. For clarity we introduce the dimensionless coupling parameter 7 = 2Vo/ftuy i . Our 
control quantity is the probability of finding the electron in the left well 



a. Localizing the electron in the left well In fig. ([T]), the time scale is given by the renormal- 
ized characteristic time T . The red dashed lines denote the optimal pulse angles obtained from 
the TLS model, which agrees well with the numerical result. In the case of full coupling fig. <[T]c), 
an initial pulse 
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FIG. 1 . The time-dependent probability Pl as a function of the applied laser pulse angle a = 2fJ, ^ p for 
different values of the rescaled coupling (a) 7 := 2Vq/ hu) p h = 0, (b) 7 = 0.6, and (c) 7 = 1. 




FIG. 2. The time-dependent probability Pl, if additionally to the localizing laser pulse with a = it/2 
a constant electric field E$ (measured in units of 2fj,/Ml) is applied, for different values of the coupling 
parameter: (a) 7 = 0.6, and (b) 7 = 1. 



with any amplitude can not totally localize the electron, but Pl still reaches values around 
90%. In order to obtain full control over the electron in this scenario as well, one can use a 
constant electric field, as explained above. We apply the unmodified laser pulse with a pulse angle 
a = 7i / 2. Fig. shows the effect on the dynamical localization. The constant electric field E 
in fig. ([2]) is rescaled in units of 2///M] so as to compare with k/VL (red dashed lines). 

b. Maintaining the localization If the electron has now been transferred predominantly to 
the left well, a periodic train of pulses with a pulse angle a = n — A/3 2 for zero coupling leads to 
a very high degree of localization. Fig. (|3| shows again the probability -Pl(£) but with the initial 
condition P^it — 0) = 1. Since the period of the pulses has been chosen as T = 80 fs, Pl shows 
a maximum for multiples of T, even in the case of a vibrational coupling. The red dashed lines 
represents the TLS results. Note that in fig. ([3]), the averaged localization is over 90 %. 
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FIG. 3. The time-dependent probability Pi as a function of the pulse angle a of the train of laser pulses 
with the T = 80 fs for values of the coupling parameter: (a) 7 = 0.6, and (b) 7 = 1. 



Conclusions We analyzed analytically and numerically the electron dynamics in a DW cou- 
pled to vibrational modes. Mapping onto a spin-model and using the polaron transformation we 
found analytical expressions for i) times and ii) light pulse parameters that are appropriate for 
switching the system from the initially delocalized to the localized target state and to maintain 
the localization. These estimates are in excellent agreement with the full-fledged two-dimensional 
numerical propagation. We considered two cases of uniform and localized vibrational coupling, 
and also when a static electric field is applied. Even in the presence of the vibrational environ- 
ment we achieve an averaged localization of 90% which makes this system a suitable candidate 
for applications in quantum information or as an ultra-fast optical switch. 
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